Continuum description of finite-size particles advected by external flows. The efl"ect of 

collisions. 
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The equation of the density field of an assembly of macroscopic particles advected by a hydro- 
dynamic flow is derived from the microscopic description of the system. This equation allows to 
recognize the role and the relative importance of the different microscopic processes implicit in the 
model: the driving of the external flow, the inertia of the particles, and the collisions among them. 
The validity of the density description is conflrmed by comparisons of numerical studies of the con- 
tinuum equation with Direct Simulation Monte Carlo (DSMC) simulations of hard disks advected 
by a chaotic flow. We show that the collisions have two competing roles: a dispersing-like effect and 
a clustering effect (even for elastic collisions). An unexpected feature is also observed in the system: 
the presence of collisions can reverse the effect of inertia, so that grains with lower inertia are more 
clusterized. 

PACS numbers: 47.52.-|-j,45.70.Mg 



Q 
U 



> 
O 

o 
m 
o 



> 

X 
S3 



I. INTRODUCTION 

The transport of inertial particles in hydrodynamic 
flows has recently attracted the attention of many re- 
searchers H. The great number of applications of this 
topic, including for e.g. ecological problems (dynamics of 
plankton populations), geophysical processes (cloud for- 
mation, transport of pollutants in the atmosphere or the 
oceans), or chemical engineering (chemical reactions in 
turbulent or chaotic flows), constitutes the main reason 
for this interest. 

In most of these works, despite of considering finite-size 
particles, the effect of the colliding processes among par- 
ticles was completely disregarded. However, in some ap- 
plications, as exemplified by rain initiation in clouds 
the role of collisions seems to be crucial. Thus, in or- 
der to take them into account, in a recent paper Q we 
have introduced a simple model of a granular material, 
N particles subject to mutual collisions (elastic or inelas- 
tic), advected by a two-dimensional chaotic flow (gravity 
is not considered). The model equivalently describes in- 
ertial particles colliding among them and immersed in a 
chaotic flow, with the density of the particles higher than 
that of the fluid. A novel result found in 3] is that col- 
lisions may strongly modify the scenario of the so-called 
preferential concentration [l| by which particles in turbu- 
lent flows, in the absence of collisions, tend to aggregate 
in specific spatial areas. 

In this Paper we derive from first principles the equa- 
tion for the macroscopic density of grains (we indistinctl 
speak of grains or colliding particles) for the model of 
(but with a flow not neccesarily chaotic). Its properties 
are discussed and, in particular, we show how it can ex- 
plain two different effects of collisions: dispersion and 
clustering. In adition, the macroscopic equation helps 
to understand another surprising feature of the model. 



observed in the Direct Simulation Monte Carlo (DSMC) 
simulations, that we call reversed clustering: in the pres- 
ence of collisions lower values of inertia can induce more 
aggregation of the particles. 

The article is organized as follows: in the next section 
we present our model and derive an explicit expression 
for the velocity of the grains. In Sec. IIIII we obtain the 
evolution equation for the density field of particles. Next, 
in Sec. II VI we present some of the most relevant features 
of this density equation; in particular, we discuss the 
relative importance of its different terms. In Sec. we 
numerically check our theoretical results, and finally, in 
Sec. IVII we summarize the Paper. 



II. EQUATIONS OF MOTION AND 
STOCHASTIC TREATMENT OF COLLISIONS 

The equations of motion of N granular particles driven 
by an external velocity field u(r, t) are given by 

^ = -4-(v,(t)-u(r,(t),t)) 

+7Ef=i [KW - V, W) • hM i^^Amt - i^)(l) 

v.W, (2) 



dv^{t) _ , 

dt 



where i = 1, ...,iV, is the velocity of particle i, its 
position, fp is the Stokes time, 7 = (1 -I- r)/2, r e [0, 1] 
is the restitution coefficient (r < 1, with r = 1 for elastic 
collisions), iiij{t) is the unitary vector joining the cen- 
ters of particles i and j at time t (in the following we will 
also use the notation n(r, r') to denote a unitary vector 
joining vectors r and r'), 6 is the Dirac delta, and with 
t'lj we denote the times at which particles i and j make 
their fc-th collision. The first term in eq. indicates 
the inertia of the particles, and the second one shows 
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how the velocity of a particle is modified because of the 
collisions with the rest. Without collisions eqs. (|1I2|I are 
the equations of motion of N spherical particles where 
the Bernouilli term, the added-mass term, the Basset- 
Boussincsq history integral term and the Faxen correc- 
tions are neglected. Thus, the advection of the particles 
is modelized taking only the Stokes drag as the relevant 
term in the Maxey and Riley equations U . This is a con- 
sistent approximation when the particles of the granular 
system are much heavier than those of the fluid , which 
is the case that we are considering. It is also fundamental 
to mention that the relevant time scales of the problem 
were identified in Q to be the typical time of the flow, 
T, the Stokes time, Tp, the mean collision time, Tc, and, 
in the case of a chaotic flow, the time rj, given by the 
inverse of the Lyapunov exponent of the flow. 

In order to obtain a close expression for Vi(t) wc first 
non-dimensionalize time and velocities taking as unit 
scales the typical time T and typical length L of the flow, 
i.e. t — > t/T, u — > uT / L and v — > vT/L. Then Vi(i) will 
be obtained as a formal expansion in the non-dimensional 
parameter Tp — From eq. ^ on has straightforwardly 



collides. 



v,(t) = u(r,(<),i) 



'^P dt 



(3) 



then this expression of is itself substituted in the r.h.s 
of lO and one obtains to order 



Vj(t) = u(ri(i),t) - Tp 



i k 

2^ /+^ , o/,3 



+ T^mit) + e(Tp^). 

If one writes down explicitly rji (t) 



(4) 



+ 7 EE(d^Hr^)-"(^^-))'^(^ 



+ 5(t-4)^(u(r,)-u(r,))) 



E E - 4)(E E(^' - ^')'^(* - 



j k 



I m 



EE(^.--^")'^(^-^'«))' 



(5) 



n p 



one can see that it is a very complicated expression tak- 
ing into account collisions among different particles. Ne- 
glecting the first term in eq. Q one can assume that it is 
some kind of effective stochastic term acting at the col- 
liding times (this hypothesis will be checked in Sec. IV)l . 
i.e., it denotes a random kick on a particle every time it 
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j=i k 



(6) 



where Cj(t) is a Gaussian white noise with zero mean 
and {Cjit)Crit')) = 2D'5{t ^ t')5jk5ira (the superindex 
denotes vector component). 

It is now easy to obtain that {r)i{t)) — 0, and the 
correlations 

= 2DS,,SkiS{t-t')n{r,it),t), (7) 

where n{ri{t),t) is such that n{r,t)drdt is the average 
number of collisions that a particle positioned in r un- 
dergoes at time t. D' and D" are constants absorbed at 
the end in the definition of D, which is another constant 
of order a'^ jT^ (tr is the diameter of the particles). This 
last statement comes from comparing the full expression 
at order in © with Q. 

Finally, denoting — yj D n{ri{t),t)r]i{t) we arrive 

to 



Vi{t) = u{r i{t),t) - Tp 



du 

'dt 



i k 

+ Tl^Dn{v,{t),t)nm + Q{Tl), (8) 

where rj is a Gaussian white noise with zero mean and 
correlations {vHtWjit')) = 25ki5ij5(t - t'). 

At this point we make the hypothesis that eq. (|HJl is 
valid for any r^; in some sense, the higher order terms, 
<d{Tp), just renormalizc the diffusion coefficient. This is 
supported by the fact that for large Tp the noise term 
in eq. (jSJ dominates, which is consistent with the fact 
that in eq. for fp large the dynamics of the particles 
is mainly driven by collisions. Thus, this hypotheses fi- 
nally states that the net effect of the colliding processes 
is properly modelized by the noise term appearing in 
This will also be checked in Section V. 



III. EVOLUTION EQUATION OF THE 
DENSITY FIELD 



Our aim is to obtain the evolution equation of the den- 
sity field of particles. For this we closely follow 0,0 and 
define p(r, t) 
Calculus |7|] 



Y!^=i Pi (i", i) = {t) - r) • Using Ito 



^^ = ~V.(p.u(r,i)) + r,V.(p.f ) 

+r4f V2(n(r, i)p,) - 7rpV • [p, Ef=i E. 
[(u(r,t) - n{vj{t),t)) ■ n(r,r,)]n(r,r,) 5{t - t%)] 



(9) 
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We assume now that the local free time between colli- 
sions Tc{r,t) is everywhere smaller than any other char- 
acteristic time of the system, and that the mean free 
path Xc{r,t) is also smaller than any other characteristic 
spatial scale. With these hypothesis we will proceed by 
integrating every term in jnj over small space-time cells 
AVAt (with this integral normalised by this spatiotem- 
poral volume) such that rc(r, i) < At < min{rp, 1, r^}, 
and A^(r, i) < AV < A^, with A the minimum typical 
space scale of the system. Assuming that the fields p and 
u are constant over the above space-time cells, the inte- 
gration over AVAt (normalised with this volume) and 
the summation over index i of the l.h.s in @, and of the 
first two terms in the r.h.s can be easily calculated. 

More complicated are the other terms. Let us study 
them in detail beginning with the diffusion term (the 
third in the r.h.s of J^): 

AT^r^*d^'r^''rfr'I^V^(p.(r',^') n(r',t')) ^ 

ri D 



2 Al^At 



E,V2(p,(r,i)/dt'dr'n(r',0) 



T„ D 



^V2p(r,t)2 



(10) 



The first equality arises from the assumption that pi is 
constant in AVAt, and in the last one we have used that 
/ dt'dr'n{r\t') w p(r,t)^^, which follows from the 
above mentioned definition of the quantity n{r, t)drdt. 
Note also that we have assumed that Tc is constant but 
all our results can be extended to a Tc dependent on the 
coordinates and time. 

Next we proceed with the collision term, the fourth 
one in the r.h.s of Eq (|5J|. Using the notation 
s(r',r",t') = [(u(r',t')-u(r",i'))-n(r',r")]n(r',r"), 
and, again, that pi is approximately constant, the sum 
over i and integration over AVAt can be written 



Then we approximate I]<i;> s(r',rj(tj^),t*j) « 
^s{r' ,rj{t),t) because ^ is approximately the number 
of colliding events in the time interval At, and the sum 
over < j >, with the help of the properties of the Dirac 
delta function, can be written as 

E<j>s(r',r,(i),i) = 
J2j j^Hhi I"' + t)6ir, (t) - (r' + an))<7, 

(13) 

so that now we do not have to manage with a restricted 
sum. Moreover, the sum of the delta's is the definition 
of p, and we have that H13|) is just 

cr /|-|^^ dn s(r',r' -I- a-n,t)p(r' + crn) 
= a2 /|_|^^ dn[{n ■ Vr'u(r', t) ■ n]np(r', t) + e{a^); 

(14) 

the equality comes from considering a Taylor expansion 
in powers of a. Therefore, the collision term takes the 
form (assuming again constancy of the density and ve- 
locity fields in AV^) 



-Vr 



Pir,t)' 



'1=1 



(in[(n • Vru(r, t) ■ n]n 



_(15) 

Finally, taking the average over the noise as in ^ , the 
last term in ^ disappears, and noting p{r,t) = {p{r,t 
we obtain the evolution equation for this density field 

-V.(pu(r,t))-KrpV.(pf ) 



dt 



r„ D 



• [P^ dn[(n • Vru(r, t) ■ n]n] . (16) 



E,,feS(r',r,(i'),mt'-4)]«-E 



p.(r, t) dt' j^y dr' E,- s(r', r,{t'), t')5{t' - t%)lll) 
Then, evaluating the time integral we obtain 



^ AVAt 



Vr-[A(r,i) / dr'^ ^s(r',r,(4),t 



<j> <k> 



(12) 

Here the X]<fc> uidicates a restriction in the sum to the 
colliding times t^j in the time interval [t, t + At] , and the 
notation X]<j> restricts the sum to the particles whose 
position is given by rj{t'^j ) = r' -I- aHj, where, as already 
indicated, a is the diameter of the particles and Uj is a 
unitary vector. This calculation deserves some clarifica- 
tions: the restricted sum to the colliding times appears 
because of the delta function, and the sum restricted to 
the particles with the position as mentioned comes from 
the fact that, right at the collision event, two particles 
are separated by a distance crfi, where n is an arbitrary 
unit vector. 



Here we have used the approximation (p^) « (p)^, since 
we expect that the densities in a point of space and time 
are likely to be uncorrelated for different realizations of 
the noise. 

It is very important to clarify the meaning of the noise 
average performed to obtain the p field. This has been 
done following the arguments of j]^. In our case, taken 
the average over the noise esentially means smoothing out 
the higher order corrections to the velocity of the par- 
ticles coming from complicated collision processes that 
take place in the time interval At. 



IV. RELATIVE IMPORTANCE OF THE TERMS 
IN THE DENSITY EQUATION, AND REVERSED 
CLUSTERING 

Some relevant features of Eq. (|16|l come inmediately 
to light. First, the inelastic character (the value of 7) 
of the collisions is almost irrelevant at this level of de- 
scription. The presence of the external driving fiow 
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turns negligible the difference between elastic and in- 
elastic colliding processes. Also, the collision-induced 
diffusivity is, as expected, dependent of the density of 
particles. This is clearly recognized by rewriting it as 

:|-^V2(p2) = lL£v(pVp) = V(L>(p)V/9), where we 

D 

have defined D{p) = 

Most of the results of the discrete model (|1I2|I that 
were presented in Q can be understood in the frame- 
work of Eq. p6(l . Let us put a label on each term of the 
r.h.s of lO: the first is the pure advective term, which 
we label with the letter a, the second is the inertial term, 
b (of order r^), which is responsible for preferential con- 
centration in the absence of collisions [l| , the third is a 
collision-induced diffusive term, c, of order t^D/tc^ and 
the last is a collision-induced term, d, of order Tpu"^ /tc, 
that could also induce clustering. 

In the absence of collisions (tc — > cxd) the last two terms 
vanish and there is only competition between the advec- 
tion, a, and the clusterizing inertial term b: for e.g. when 
Tp >~ 1 it appears that inertia dominates and the sys- 
tem is strongly clusterized. On the other side, in the limit 
Tp ^ (6 vanishes), it is the d term that can induce clus- 
tering (when Tc Tp); however the term c may eventu- 
ally become stronger than d giving rise to a homogenous 
distribution of particles. The unexpected feature emerg- 
ing from this analysis is that the diffusive homogenizing 
term c dominates at high values of Tp, thus reversing the 
effect of inertia that, in collision-less systems, enhances 
concentration. In the presence of collisions, therefore, we 
have a reversed clustering phenomenon: for small values 
of Tp, i.e. inertia, the density is more clusterized than at 
larger values. 



V. NUMERICAL SIMULATIONS 

We have checked numerically our system in two spatial 
dimensions by means of DSMC simulations of the hard 
disks model (similar to those performed in 3]) and by 
a numerical study (a Lax integration scheme) of eq. . 
The DSMC is a well established algorithm 9] that allows 
to simulate gaseous systems in a rapid and efficient way, 
with the assumption of negligible correlations at short 
range. We have used a variant of DMSC that takes into 
account the Enskog factor due to high density correc- 
tions, in order to obtain more accurate simulations of the 
clusterized situations. In Appendix we give details about 
this simulation scheme. First of all, we have verified the 
stochastic approximation done in eq. ©, using a simple 
velocity field given by Ux{ri{t),t) = U cos{2TTyi{t)/ L), 
Uy — 0, with Yi(t) = {xi{t),yi{t)). We calculate the 
quantity given by the l.h.s of |(SJ) minus the first three 
terms in the r.h.s, cumulated for a little time period At 
(greater than the simulation time step, but shorter than 
the mean free time) and divided by Tpy/n. We have per- 
formed this for both cases: large and small Tp, in order 
to check also the discussion below eq. One can ap- 
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FIG. 1; Distributions of the residual term (left: x compo- 
nent, right: y component) of eq. ||HJ which is assumed to be 
a white Gaussian noise. The form of the distribution and 
the self-correlation (displayed in the insets) confirm this as- 
sumption. The flow used in this simulation is a sinusoidal 
shear, i.e.: Ux{ri{t),t) — U cos{2Tiyi(t) / L), Uy = 0, where 
ri(t) — {xi{t),yi{t)). For the upper plot (large Tp) the time 
step used in the DSMC is At = 0.0001, size of the DSMC 
cells is 0.6 X 0.6, mean free time is 0.01, mean free path 5, 
length of the time period used to calculate r; is At = 0.005, 
Tp = 1, r = 1, A = = 1000, a = 0.3, U = 10000 so that 
Tp = ^ 300. For the bottom plot (small Tp) At = 0.01, 
DSMC cells are 2x2, mean free time is 0.6, mean free path 
18, At = 0.5, Tp = 1, r = 1, A = = 10000, a = 0.6, 
f/ = 50 so that Tp = 0.5. 



preciate in figure ^ that in the two cases this quantity is 
a Gaussian white noise. 

Typical patterns for the distribution of particles ob- 
tained with DSMC at large time can be seen in Fig. El 
The flow we have used here is the cellular flow derived 
from the streamfunction PI ^p{x, y,t)^U sin(27r/L(a; + 
Bq cos{ujt))) sm[2'K / Ly) , with Bq and uj the amplitude 
and frequency of the temporal perturbation, respectively. 
In this flgure we just show a small part, around 2% of 
the entire system, in order to better appreciate the clus- 
tering areas. For the sake of clarity we also show, in 
red lines, the streamlines of the flow for Bq — Q. We 
have chosen two opposite situations. In figures |2K and 
|2J3 we have studied a strong inertia case (jp > 1) with 
the panel [2}3 corresponding to the same case as|2K but 
without collisions. Note that in the absence of collisions 
b > a and preferential concentration can be recognized 
(panel [5)3), while with (elastic) collisions we have that 
c^d^b^a so that diffusion dominates and the dis- 
tribution of particles in[2K is homogeneous. Thus, in this 
case we observe the dispersing effect of the collisions. In 
figures I2P and[5j3, instead, we have considered a weak 
inertia case, Tp < 1, where |2f) is the case with no colli- 
sions. Now, without collisions we have that a > b and 
the flow homogenizes the distribution of grains (see|2f)). 
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With collisions, instead, we have d > a > b > c and colli- 
sion induced clustering is observed Therefore, this 
case shows the clustering effect of the collisions. 




the Pm {n) resembles a Dirac Delta distribution centered 
around this constant value. 

The results shown in figs. Inland |31 clearly indicate the 
relevant role of the collisions for inertial particles im- 
mersed in a flow, verifying as well the important and 
unexpected phenomenon of reversed clustering that ap- 
pears on the basis of the inspection of terms c and d of 
equation (|16ll . Briefly, in the absence of collisions parti- 
cles with inertia tend to aggregate, and the largest is the 
value of Tp the more compact is the aggregation of par- 
ticles Q. This can be seen comparing figs. |33 and|2p). 
However, when collisions are taken into account this last 
statement can be wrong, and just the contrary occurs: in- 
ertia is decreasing but aggregation increases (again this is 
what happens if we compare figs.|2K and|2P). Note that 
in these last two cases we have that Tc < Tp so that col- 
lisions dominate the dynamics. Equation perfectly 
reflects this situation: the main difference between the 
two cases is originated by the value of Tp that influences 
the relative importance of terms c and d, as well as the 
relative importance between a and b in the absence of 
collisions. 

It is important to note that this effect may also appear 
in segregation proccesses, i.e., in systems with particles 
of different sizes or densities. Obviously, these particles 
have different values of Tp and the same scenario just 
commented emerges. Finally we put in evidence that the 
inelasticity of collisions (typical of granular materials) 
here just plays the role of slightly enhance spatial corre- 
lations, i.e. clustering. This is not taken into account by 
equation (fTC|l but can be observed in the DSMC simula- 
tions. 



FIG. 2: Instantaneous density patterns, obtained with 
DSMC, for the 2d cellular flow; N = 500000 and L = VN. 
B and D are without collisions. Frames A, and B: U — 30, 
Bo = 1 = 1 (r/ ~ 0.4), Tp = 10, o- = 0.3, r = 1, Tc ~ 0.02 
(only for A). Figures C and D: [/ = 110, Bo = 1, w = 10 
(r/ ~ 0.3), Tp = 0.003, a = 0.5, r = 1, Tc ~ 0.0002 (for C). 
come from DSMC. on Note that a small part (~ 2%) of the 
whole system is portraid. Red lines represent streamlines of 
the flux with Bo — 0. The inset in B shows these for the 
entire spatial domain. 

The numerical comparison of the DSMC with the con- 
tinuum equation is shown in Fig. |31 Here we measure 
the PM(n) function, which gives the histogram (normal- 
ized to unity) of the number of boxes containing n parti- 
cles after dividing the system in M boxes. Note that as 
the clustering is stronger the deviations of Pm (n) from a 
Poissonian (homogenous distribution) are more evident. 
On the right wc plot the function for the DSMC patterns 
of flgs. and |33 (upper right panel) and flgs. |2P and 
(lower right). On the left, we show the Pm("-) for 
the distribution calculated with the continuum equation 
and the same value of the parameters. Note that now 
a homogenous distribution of particles corresponds to a 
constant density and that is the reason why in this case 



Continuum Description Microscopic Description 




FIG. 3: PM(n) function for the distribution of particles ob- 
tained from DSMC (right panels) and for the density of par- 
ticles from the continuum equation (1161 (left) for the same 
parameter values. Upper right are the cases in |2K ^nd |2J3, 
that is, with high inertia, and the lower right corresponds to 
the|^ and|5f) panels (low inertia). 
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VI. SUMMARY 

The evolution of the density field of a large number of 
colliding finite-size particles driven by an external flow 
has been derived. It is important to note that no premise 
has been perfomed on the flow so our results are equally 
valid for laminar, chaotic or turbulent flows. However, in 
the numerics we have used a chaotic flow just for its inter- 
est and simplecity, and also to compare with the results 
of our previous work . In this equation one can recog- 
nize two, in principle, competing roles of the collisions: 
clustering and dispersion. The relative strengths of the 
important time scales involved in the system will flnally 
determine if the particles tend to cluster or, on the con- 
trary, are uniformly distributed in the space. Collisions 



may even invert the typical scenario found for inertial 
particles, showing the reversed clustering phenomenon. 
We have also presented a numerical study of the con- 
tinuum equation and compared with DSMC simulations 
of the system of discrete particles. The good agreement 
gives support to our theory. 



VII. ACKNOWLEDGMENTS 

We acknowledge fruitful discussions with Umberto 
Marini Bettolo Marconi. C.L. acknowledges support 
from MCyT of Spain under project REN2001-0802-C02- 
01/MAR (IMAGEN). 



[1] T. Elperin, N. Kleeorin, and I. Rogachevskii, Phys. Rev. 
Lett. 77, 5373 (1996); E. Balkovsky, G. Falkovich, and A. 
Fouxon, Phys. Rev. Lett. 86, 2790 (2001); T. Nishikawa, 
Z. Toroczkai, C. Grebogi, and T. Tel, Phys. Rev. E 65, 
026216 (2001); C. Lopez, Phys. Rev. E 66, 027202 (2002); 
I.J. Benczik, Z. Toroczkai, and T. Tel, Phys. Rev. E 67, 
036303 (2003). 

[2] G. Falkovich, A. Fouxon and M.G. Stepanov, Nature 

419, 151 (2002). 
[3] C. Lopez and A. Puglisi, Phys. Rev. E 67, 041302 (2003). 
[4] M.R. Maxey and J.J. Riley, Phys. Fluids 26, 883 (1983). 
[5] M.R. Maxey, J. Fluid Mech. 174, 441 (1987). 
[6] U.M.B. Marconi and P. Tarazona, J. Chem. Phys. 110, 

8032 (1999). 

[7] D.S. Dean, J. Phys. A: Math. Gen. 29, L613 (1996). 
[8] In two spatial dimensions the integral in the col- 
lision term can be performed and one obtains 

[9] G.A. Bird, Molecular gas dynamics (Clarendon Press, 
Oxford, 1976). 

[10] T.H. Solomon and J.P.Gollub, Phys. Rev. A 38, 6280 
(1988). 

[11] C. Cercignani, R. lUner, and P. Pulvirenti, The Math- 
ematical Theory of Dilute Gases (Springer-Verlag, New 
York, 1994). 

[12] W. F. Carnahan and K. E. Starhng, J. Chem. Phys. 51, 
635 (1969). 



Appendix: the DSMC 

Direct Simulation Monte Carlo, also known as Bird 
Method 01, is a simulation scheme used to solve Boltz- 
mann equations. With some conditions and in well de- 
fined limits its results are proved to converge to the 
solution of the Boltzmann equation for a gas of hard 
spheres jllj . The algorithm consists in two main steps: 
a) the free flow, and b) the collisions. Every time step of 
length At (lesser than any characteristic physical time. 



e.g. the collision time) all the particles are moved freely, 
i.e. disregarding possible collisions and then the collision 
procedure is applied: the system is divided into cells of 
linear size lesser than any characteristic physical length 
(e.g. the mean free path) and particles in the same cell 
may collide randomly. For each cell c a number of 
couples of particles is chosen with probability propor- 
tional to their scattering section, i.e. in this case pro- 
portional to their relative velocities. The number nic is 
calculated as rric = tOcAt, where cOc is the average col- 
lision frequency estimated in the cell c, by assuming a 
Gaussian distribution of velocities with variance given by 
the actual variance of velocities of the particles inside the 
cell. As it is known, the Boltzmann equation is a correct 
description of a gas of hard particles only in the dilute 
limit TV ^ oo, (7 — > with flnite Na'^~^ {d the space 
dimension). When the gas cannot be considered dilute, 
correlations arise in the form of an enhancement of the 
collision frequency and in complicated excluded volume 
effects. It is accepted that at not too high packing frac- 
tion, the so-called Enskog correction to the Boltzmann 
equation gives a sufficient description of these effects. In 
the Boltzmann equation this correction appears as a sim- 
ple multiplicative term in front of the collision integral, 
which is equivalent to an increase of the collision fre- 
quency. The Enskog correction is usually taken to be the 
static correlation function g(r) evaluated at contact, i.e. 
r — a, for which approximated forms (dependent upon 
the local volume fraction (p = Nc'K{cr /2Y /Vc, with Vc the 
area of a cell and N^. the number of particles in the cell) 
are available. We (in d — 2) have used the following 
form [l3: 
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